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Pulsar winds are the ideal environment for the study of non-linear electromagnetic waves. 
It is generally thought that a pulsar launches a striped wind, a magnetohydrodynamic 
entropy wave, where plasma sheets carried along with the flow separate regions of alter- 
nating magnetic field. But when the density drops below a critical value, or equivalently 
for distances from the pulsar greater than a critical radius, a strong superluminal wave 
can also propagate. In this contribution we discuss the conversion of the equatorial 
striped wind into a linearly polarized superluminal wave, and we argue that this mode 
is important for the conversion of Poynting flux to kinetic energy flux before the outflow 
reaches the termination shock. 
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flows. 
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1. Introduction 

Pulsar winds are relativistic outflows, consisting mainly of electrons and positrons. 
They are thought to be Poynting-flux dominated at launch, as quantified by the 
magnetization parameter a which is the ratio of Poynting to kinetic energy flux, 
with cr S> 1. Far away from the light cylinder, which is at radius tlc — c/uj with 
w the pulsar's rotational frequency, the magnetic field is predominantly toroidal 
with an amplitude falling as 1 /r and reverses sign across a corrugated current sheet 
that separates the two magnetic hemispheres. The wind propagates radially as a 
transverse subluminal wave with speed equal to the bulk speed of the outflow, or 
equivalently Lorentz factor F. The opening angle around the equator of the region 
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of alternating magnetic field depends on the misalignment of the magnetic and 
rotational axes of the pulsar, denoted by a. 

In the absence of dissipation, the wave propagates outwards from the light cylin- 
der while remaining Poynting-dominated"'^. However, after crossing the termination 

o 

shock, the outflow's magnetization is low - . Where and when conversion of Poynt- 
ing flux to particle energy happens is still not clear. Magnetic reconnection in the 
current sheets separating regions of opposite magnetic field polarity in the wind is 
a possible dissipation mechanism, however this process is much too slow to account 
for the conversion in isolated pulsars like the Crab, unless the pair injection by 
the pulsar is much larger than generally assumed^^. Another possibility is driven 
reconnection at the termination shock, which has been demonstrated to dissipate 
the striped wind's alternating field efficiently in some cases. This, however, requires 
even higher particle injection rates than reconnection in the current sheets^™. 

A solution can be provided in the form of non-linear superluminal waves, i.e. 
electromagnetic waves of large amplitude with phase speed > c. These can 
propagate outside of a critical radius r^, at regions of lower particle densitjEl (as 
a contrast to the magnetohydrodynamic wind which needs a minimum particle 
density) . The conversion of the striped wind to a superluminal mode is possible if 
Tcr is smaller than the shock radius. 

In the following we will investigate non-linear superluminal waves of linear po- 
larization, which are relevant in a region of opening angle a around the equator of 
a pulsar wind. We identify the superluminal modes into which a striped wind can 
convert by imposing "jump conditions" across the (thin) region where conversion 
happens. We will assume throughout that the waves can be considered as locally 
plane and purely transverse. These are excellent approximations far away from the 
light cylinder (at p — r/ri^c ^ 1). The radial direction, which is the direction of 
propagation of the wave, can be identified with the x— direction in a cartesian system 
of coordinates. The surviving field components are then Ey and B^, corresponding 
to the azimuthal magnetic field and the polar electric field. 



2. Transverse superluminal waves: the two-fluid approach 

The equations governing the evolution of a cold, two-component plasma are the con- 
tinuity equation and the equations of motion for each species, along with Maxwell's 
equations. We are looking for periodic solutions in a special frame of reference 
that propagates with velocity (? jv^ with respect to the laboratory frame (in the 
laboratory frame, which coincides with the pulsar frame, the wave phase speed is 
i;^). We will call this special frame "homogeneous "or "H-frame". The advantage 
of this formulation is that in the H-frame all space dependence vanishes and only 
time dependence remains, and the phase is just woi with wq the wave frequency!^. 
From now on quantities in the H-frame will be unprimed, while quantities in the 
laboratory frame will be primed. 

In the H-frame, the conservation of magnetic flux (V • B = 0) is automatic, and 
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it can be seen immediately from Faraday's law that the magnetic field is constant. 
From Coulomb's law and the continuity equation we get n+7+ = n_7_, where a 
plus(minus) subscript denotes positrons(electrons). Furthermore, the Lorentz fac- 
tors are equal 7+ = 7_ and the other components of the dimensionless four- velocity 
of the two species are related by: 



Because of these equations we can drop the subscripts and solve the system in 
terms of the positron fluid quantities. It can be showrM^ that all quantities can be 



Ux+ = Ux~ Uy+ = -Uy- Uz+ = -Uz- ■ (l) 

ri pts 
.rPl 

expressed in terms of the electric field normalized to its largest amplitude y = E/Eq: 

7-7o + 27o(l-2/')/g (2) 

= y^7o'-l + 47oA(l-y)/g (3) 

Uy = ±^^^-ul^l (4) 

Uz^O. (5) 

In the above, quantities with a zero subscript express the initial condition at zero 
phase, and the field reaches its largest amplitude at phase zero, yo — ^- Also, 
X = B/Eo is the ratio of the (constant) magnetic field to the electric field amplitude, 
and q is double the ratio at phase zero of the energy density in the particles to that 
in the electric field: 

1 = ■ (6) 

Eq 

The electric field y can be computed from a first-order differential equation, which 
in the general case (A ^ 0) has no solution in closed form and has to be numerically 
integrated. The phase-averaged values of the above quantities, then, can also be 
calculated numerically (for details see Ref . IT^ . 



3. Conserved quantities 

There are four conserved quantities in a pulsar wind, in the plane wave approxima- 
tion. These are the phase-averaged x— components of the particle, energy, momen- 
tum and magnetic flux densities. We will denote the particle flux by (J'), while the 
energy and momentum flux densities are the ^T'"^^ and {T'^^'^ components of the 
combined stress-energy four-tensor of fields and particles. Finally, the magnetic fiux 
density can be replaced (from Faraday's law) by the phase-averaged electric field, 
so that the last conserved quantity can be written as (£") /{iir). 

From these four quantitites we can construct three dimensionless ones, defined 

as 
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Fig. 1 . On the left the index of refraction is plotted as a function of dimensionless radius _R, for 
four different cases of b. On the right, the phase-averaged Lorentz factor of the particles is plotted. 
We see that for low b, {7') ~ /i, which means that almost the whole energy of the wind is now 
carried by the particle component of the flow. 

V = ^ , ' (8) 

V=^^ (9) 
' 47rmc (J') ■ ^ ' 

The first of these, /i, is the energy per unit mass that each particle would have, 
if the entire energy flux were carried by the particles. The three quantities (I?])-® 
are conserved across a layer that marks the transition between the striped wind 
and the superluminal wave, and therefore they have to be expressed in terms of the 
quantities characterizing these modes. For the striped wind we have: 

fi = r{i + <j) (10) 

= [r2(l + a) - (1 + a/2)] (F^ - l)-i/2 (n) 
jy = b^Ta^ (12) 

where a is the wind magnetization and b is the normalized DC component of the 
magnetic field. This quantity is zero at the equatorial plane, where regions of al- 
ternating field are equal in width and cancel out in the phase average, and rises 
to its maximum value, 5 = 1, at latitude equal to the inclination angle a, where 
the alternating field component disappears. The same quantities /x, v and rj can be 
calculated for the superluminal wave. Since our treatment of these modes took place 
in the H-frame, we calculate the quantities (J), {T°°), (r"i), (T") and {Ef /{An) 
and transform them to the laboratory frame using the velocity /v,f,, as explained 
above. 

Equating Eqns. ©-(jl]) calculated for the superluminal mode to Eqns. (fTO|)- 
(fT2|) . one can find "jump conditions" for the transition from the striped wind to the 
superluminal mode, for certain values of F, a and latitude, or equivalently b. We 
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Fig. 2. The critical surfaces for three values of the inclination angle a, along with two example 
shock surface, normalized to equatorial radius peq.ts two and five times larger than the equatorial 
critical radius Peq,cr. The pulsar is at the point (0, 0) and one half of the poloidal plane is shown. 
The variable r is normalized to the equatorial value of R. 

can then re-introduce the conservation of particle flux into our system of equations, 
in order to find the radial dependence of these jump conditions^. 

4. Results and discussion 

In Fig. [T] we have plotted as a function of normalized radius R — pfx/a^ the inverse 
hyperbolic tangent of the refractive index c/v^ of the wave into which a striped 
wind with F = 100 and a — 100 can convert. The quantity is a dimensionless 
parameter given by: 



"^ = V;^rff7 = "^°'^"^' 

where by dL/dCl the distribution of the pulsar's luminosity in solid angle is meant 
and 6 is the colatitude (i.e. the polar angle). We have assumed that the distribution 
of the pulsar's luminosity follows the prescription 

f^Losin^. (14) 

and separated the angular dependence from al. 

As we can see in Fig. [1] there is a minimum normalized distance i?cr outside 
which superluminal waves can propagate. For 6 = 0, Ra — 1, but it rises as one goes 
to higher latitudes. On the plot on the right of Fig. [T] the phase averaged Lorentz 
factor of the particles, (7') is shown for the same values of b. From this plot it is seen 
that during the conversion from one wave to another, energy is transferred from 
fields to particles. In the striped wind the Lorentz factor of each particle is equal to 
the bulk Lorentz factor of the outflow, which in the plotted example is F = 100 <C 
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However, in the superluminal wave, for low values of 6, (7') ~ /z, which means that 
almost all the energy of the outflow is now carried by particles. Even for values of 
b close to unity, (7') F, meaning that the particles get accelerated in the whole 
wind, and the ratio of the Poynting to the kinetic energy flux decreases during the 
transition from the subluminal to the superluminal regime. 

In Fig. [2l the geometrical critical surface is shown (radius normalized as p ^■ 
p/i/flLo), outside which superluminal waves can propagate, for different values of 
the inclination a. Superimposed on the critical surfaces is the example of a shock 
surface, the shape of which is taken from Ref. [131 The normalization of the shock 
radius in comparison to the critical radius in Fig.[2]is arbitrary. Generally it depends 
on the individual pulsar in question. For isolated pulsars like the Crab and the 
Vela, the shock (located at ^ lO^rLC and ~ lO^rLc respectively^^) is much farther 
out than the critical surface (located at lO^rLc and ~ lO^rLC respectively) for 
practically the whole of the wind. For pulsars in binary systems, though, where the 
shock is pushed closer to the pulsar by the outflow from a companion stai^^, it 
is possible that the critical surface is partly or entirely outside the shock, so that 
no superluminal waves can propagate. In highly eccentric systems, the possibility 
arises that the termination shock migrates during one orbit from a location inside 
to one outside the critical surface, so that the characteristics of the wind just before 
the termination shock are orbitally modulated during one period of the system. 

We conclude that the conversion of the striped wind to a superluminal mode is 
possible in many objects, especially in isolated pulsars, and it is accompanied by an 
acceleration of the particles of the outflow to (7') ^ F, which might influence the 
physics of particle acceleration and the structure of the termination shock. 



References 

1. J. G. Kirk and I. Mochol, Ap. J. 729, 104 (2011). 

2. M. J. Rees and J. E. Gunn, Mon. Not. R. Astron. Soc. 167, 1 (1974). 

3. F. V. Coroniti, Astron. Astroph. 349, 538 (1990). 

4. Y. Lyubarsky and J. G. Kirk, Ap. J. 547, 437 (2001). 

5. J. G. Kirk and O. Skjaeraasen, Ap. J. 591, 366 (2003). 

6. Y. Lyubarsky, Mon. Not. R. Astron. Soc. 345, 153 (2003). 

7. Y. Lyubarsky and J. Petri, Astron. Astroph. 473, 683 (2007). 

8. L. Sironi and A. Spitkovsky, Ap. J741, 39 (2011). 

9. A. Melatos and D. B. MeIrose,Mon. Not. R. Astron. 500.279, 1168 (1996). 

10. P. C. Clemmow, J. PI. Ph. 12, 297 (1974). 

11. C. F. Kennel and R. Pellat, J. PL Ph.l5, 335 (1976). 

12. I. Arka and J. G. Kirk, ^arXiv:! 109. 27S6v l (accepted for publication in Ap. J.) (2011). 

13. Y. Lyubarsky, Mon. Not. R. Astron. Soc. 329, L34 (2002). 

14. O. Kargaltsev and G. G. Pavlov eprint arXiv: 0801.2602v2, (2009). 

15. Chaterjee. S. et al, Ap. J. x670, 1301 (2007). 



